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Abstract We consider a late closed universe of which 
scale factor is a power function of time using ob- 
servational data from combined WMAP5+BA0+SNIa 
dataset and WMAP5 dataset. The WMAP5 data give 
power-law exponent, a = 1.01 agreeing with the pre- 
vious study of H(z) data while combined data gives 
a = 0.985. Considering a scalar field dark energy and 
dust fluid evolving in the power-law universe, we find 
field potential, field solution and equation of state pa- 
rameters. Decaying from dark matter into dark energy 
is allowed in addition to the non-interaction case. Time 
scale characterizing domination of the kinematic expan- 
sion terms over the dust and curvature terms in the 
scalar field potential are found to be approximately 5.3 
to 5.5 Gyr. The interaction affects in slightly lowering 
the height of scalar potential and slightly shifting po- 
tential curves rightwards to later time. Mass potential 
function of the interacting Lagrangian term is found to 
be exponentially decay function. 



1 Introduction 

The presence of a scalar field is motivated by ideas in 
high energy physics and quantum gravities, although it 
has not been discovered experimentally. TeV-scale ex- 
periments at LHC and Tevatron might be able to con- 
firm its existence. It is nevertheless widely accepted in 
several theoretical modeling frameworks, especially in 
contemporary cosmology of which an early-time accel- 
erating expansion, the inflation is proposed to be driven 

Burin Gumjudpai 
Kiattisak Thepsuriya 

^ThEP's CRL, NEP, Tlie Institute for Fundamental Study, Nare- 
suan University, Phitsanulok 65000, Thailand 

^Thailand Center of Excellence in Physics, Ministry of Educa- 
tion, Bangkok 10400, Thailand 



by a scalar field in order to solve horizon and flatness 
problems (Starobinsky et. al. 1980). After inflation. 



components of barotropic fluids such as radiation and 
other non-relativistic matter were produced during re- 
heating and cooling-down processes. The present uni- 
verse is also found to be in acceleration phase which 
is strongly backed up by various observations, e.g. the 
cosmic microwave background (Masi 2002[ ), large-scale 
structure surveys (Scrantan et. al. 2003) and super- 



novae (SN) type la observations (Riess et. al. 1998 



2004| |Astier]|2005[ ). The scalar field could be responsi 
ble to the present acceleration in various models of dark 
energy ( Padmanabhan et. al.|[2004 ). 

Power-law cosmology with a oc t", describes an ac- 
celeration phase for with a > 1. The model provides 
simple model of expansion as well as a convenient way 
to fit with observations. Although tightly constrained 



by nucleosynthesis (Sethi et. al. 1999 Kaplinghat et 



al. 1999), the power-law expansion with a ~ 1, could 



be considered viable at later time sufficiently long after 
matter-radiation equality era. Studies of high-redshift 



objects such as globular clusters (Kaplinghat et. al. 



T999l |Lohiy aet. al.|1997l |Dev et. al.l2002| [Sethi et. al. 



20051 , SN la data ( |Sethi et. al.|2005[]Kumar|201lD, SN 



la with H{z) data (Devet. al.|2001 Sethi et. al.[2005 



Dev et. aL]|2008 ) and of X-ray gas mass fraction mea- 



surement of galaxy clusters (Allen et. al. 2002 Zhu et. 



al. 2007) within context of power-law cosmology have 
been reported from time to time. Moreover other as- 



pects such as gravitational lensing statistics (Dev et. al 



2002 1 and angular size-redshift data of compact radio 



sources ( |Jain et. al.|[2003| of the power-law cosmology 
can be applied to constrain the model. Motivations 
of the model are various. Simplest inflationary model 
with power-law expansion can remove flatness and hori- 



zon problems leaving with simple spectrum (Lucchin 



et. al. 1985). In present universe, idea of hnear coast- 



ing cosmology (a cx t) (Kolb 1989) are proposed such 
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that it can resolve age problem in the CDM model ( Lo- 



hiya et. al.|1997 1 . On fundamental physics motivation, 
the linear coasting model can arise from non-minimally 
coupled scalar-tensor theory in which the scalar field 
couples to the curvature contributing to energy density 



that cancels out vacuum energy (Sethi et. al. 1999 



Ford|[r987| ). The model could also be a result of domi 
nation of an SU(2) cosmological instanton in the early 
universe ( [Allen |[l999l ). 

We consider a non-flat {k ^ 0) FLRW universe 
with expansion in the form of the power law function. 
The power-law expansion occurs well after matter- 
radiation equality era. Two major ingredients are scalar 
field dark energy evolving under the scalar field po- 
tential V{(j)), and barotropic (dust) fluid consisting 
of cold dark matter and baryons. We use combined 
WMAP5-HBA0-HSN d ata and the WMAP5 data alone 
( Hinshaw et. al.pOOS ) to determine scalar field poten- 
tial, values of a and other relevant parameters. Other 
related study for scalar field power-law cosmology are 



for the flat case and with phantom field (Kaeonikhom 



et. al. 2010). Previously, not within the context of 



power-law cosmology, there are a number of works in 
constructing scalar field potential directly from obser- 
vational data or in theoretical construction of potential 
under specific assumption. These are considered in var- 
ious situations, for example, construction of the poten- 



tial using SNe la dataset ( Simon et. al. 2005 1, construc- 



tion of potential assuming scaling solution (Liddle et 



al. 19991, case of scalar field negative potential with- 



out using observational data (Cardenas et. al. 2004), 



construction of scalar field potential using dark energy 



density function with non-specified expansion law ( Guo 



et. al.|2006 ) and other several works in (Huterer et. al 
1999). In this work, we also consider situation when 



there is interaction between cold dark matter and dark 



in the late FLRW universe of wCDM model with 
zero cosmological constant are considered. The time 
evolution of the barotropic fluid (cold dark matter and 
baryonic matter) is governed by the fluid equation 



since is constant, 

_ D 
a" 



(1) 



(2) 



where n = 3(1 -|- Wm) and Z? > is a proportional 
constant. For the scalar field, supposed that it is min- 
imally coupled to gravity, its Lagrangian density is 
C = (/)2/2 — V{(l>). The energy density and pressure 



(3) 



The fluid equation of the fleld describing its energy con- 
servation as the universe expands is 



■iH4)+ —V = 0. 



(4) 



Total energy density ptot and total pressure ptot of the 
mixture are simply the sums of those contributed by 
each fluid, for which the Friedmann equation in SI unit 
is 



H' = 



SttG 



"Ptot 
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It is straightforward to show that 



V 



M4c / 2 



H 
"3 



2fcc2 



6\ Dc^ 



(5) 



(6) 



energy with constant interaction rate ( |Amendola|2000[ ) . The reduced Planck mass, Mp is related to Newton's 

constant G as Mp = Hc/SttG. For dust fluid n = 3 



We see how the interaction can affect scalar field po- 
tential and we will find mass potential function of the 
interaction Lagrangian. We stress that the WMAP5 
combined dataset allows non-flat universe and dark en- 
ergy equation of state. This is because the addition of 
the SN la analysis to the WMAP5 dataset marginalize 
the equation of state parameter and spatial curvature 
at the same time. The later WMAP7 dataset is not 
suitable our study since its relevant derived parameters 
are found under flat assumption, allowing only studying 
of equation of state. 



2 Scalar-field power-law cosmology 

Here we will work with observational data in SI units. 
Two perfect fluids, the dust matter and scalar field (j) = 



{wra = 0), the last term is just —Dc^ /2a? . In power- 
law cosmology, scale factor is a power function of cosmic 
time 



a{t) = ao 



Hence the Hubble parameter is 
a[t) a 



a{t) 



t 



(7) 



(8) 



Incorporating ([7| and ([8| into the above equation, we 
obtain 



V{t) 



i2 



2 agfS"' 



(9) 
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We see that linit_j.oo V{t) = at late time. We will con- 
sider contribution of the first term (effects of acceler- 
ation rate and Hubble expansion) alone in comparison 
to total contribution when including the second (the 
curvature) and the third (dust density) terms. 



3 Cosmological parameters and observational 
data 

Using ([8| at the present time, we have a = -ffo^o- The 
sign of k depends on the sign of the density parameter 
flk = —kc^/a?H^. In our convention here, k = 1 (£7^ < 
0) for a closed universe. A: = for a flat one, and k = 
— 1 (J7fe > 0) for an open one. Allowing k to take these 
three discrete values (0,±1), present value of the scale 
factor, ao is not fixed to unity and it can be found from 
the definition of ilj. g directljQ 



ao 




(10) 



The density constant D can be found from ([2|, 

D = Pm,0 flo = ^m,0 Pc,0 Aq, (H) 

where rini,o = i^CDM,o + ^^6,o, i-e. the sum of the present 
density parameters of the barotropic fluid components. 
Pc,o is present value of critical density. Radiation and 
neutrino contributions are negligible. The values of H^, 
to, ^k,Q, i^CDM,Oi and ftb.o are taken from observational 
data. Note that although to value in the WMAP5 data 
is derived in ACDM model. This fact can be relaxed 
and we can estimably use it. This is because the ob- 
served equation of state at present is very closed to -1 
and to can not be derived if neither the equation of 
state nor a is given. Two sets of data used here are 
provided by (Hinshaw et. al. 2008). One comes solely 



from the WMAP5 data and the other is the WMAP5 
data combined with distance measurements from Type 
la supernovae (SN) and the Baryon Acoustic Oscilla- 
tions (BAO) in the distribution of galaxies. For to, Hq, 
Hb^Oj and ricDM,Oi we take their maximum likelihood 
values. The curvature density parameter ftk,o comes 
as a range with 95% confidence level on deviation from 
the simplest ACDM model. These data can be seen 
in Table ^ in which flk.o is reported as an interval of 
minimum and maximum values. We use the mean of 
each ilfc interval to represent ftk.o- Table [2] presents 
parameters derived from data in Table [l] The exponent 
a can be simply found using a — HqIq. The WMAP5 



^When considering case of dark-matcr and dark-energy interac- 
tion, we set ao = 1 for convenience. 



data can give an accelerating expansion as a > 1 while 
the combined WMAP5+BA0+SN does not. The de- 
termination of a from other studies are such as X-Ray 
gas mass fractions in galaxy clusters which favors open 
universe with a > 1 (a = 1.14 ± 0.05) (Zhu et. al 



2007 1 . Larger number of data points are needed in this 



study in order to distinct cosmic geometry in power- 
law scenario. Combined analysis from SNLS and H{z) 
data (from Germini Deep Deep Survey-|-archival data 
points) yields a = 1.31 (assuming open geometry) (Dev 



et. al. 2008). Flat and closed geometry are still vi- 



able in these studies. Compared to ( Dev et. ar]|2008 ) 
without assuming geometry, H{z) (spatial-geometry in- 
dependent) data alone gives 1.07 while WMAP5 data 
alone yields a ~ 1.01. Combined i7(z)+SN la data 



(SNLS) assuming closed geometry gives a — 2.28 (Dev 
et. al. 2008). Although Table [2] shows mean values of 
^k,o in closed geometry region, but it is clear that with 
closed geometry, our values of a do not match the a 
value from i/(z)-fSN la data. The curvature depen- 
dence in SN data analysis could give larger deviation 
in a from unity. Better SN data analysis would im- 
prove determination of spatial curvature in power-law 
cosmology. 



4 Results : non-interacting case 

Using combined WMAP5+BA0+SN dataset, the po- 
tential is 



1.03 X 10^ 

t2 



1.51 X 10^3 1.50 X 10^ 



^1.97 



t2-95 



(12) 



whereas, for WMAP5 dataset alone. 



V{t) 



1.11 X 10- 

t2 



,26 



7.67 X 10- 



,24 



4.69 X 10 



43 



(13) 



in SI units. Their plots are shown in Fig. [TJ The 
points at which the potential, its derivative, and its 
second-order derivative, are zero (tinterceptj ^maxi and 
^infloctionj respectively) are also determined, for both 
Hj. and each end of the ^k,o interval. The results are 
summarized in Table [2] After tinflcction, the potential 
from each data behaves like its first term, i.e. decreas- 
ing in its value while increasing in its slope (being less 
and less negative). The other terms quickly become 
weaker. This can be seen in Fig. [l] Since the first term 
is contributed only by H{t) (and its time derivative), 
it is dominant in the post-inflection phase. In fact, the 
convergence to zero of the potential is slower than its 



flrst term alone (see ( 12 ) and ( 13 )), because the sum of 
the last two terms consequently becomes positive before 
converging to zero. This means that the plots of each 
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Table 1 Observational derived parameters from ( Hinshawet. al.|2008 l 
Parameters WMAP5+BAO+SN WMAP5 



to 
Ho 

^b.O 

^k,0 
Pc,0 
Pb,0 
Pmfi 



13.72 Gyr 
70.2 km/s/Mpc 
0.0459 
0.231 

-0.0179 < ilkfl < 0.0081 
9.26 X 10^27 kg/m^ 
4.25 X 10^28 kg/m^ 
2.56 X 10^27 kg/m^ 



13.69 Gyr 
72.4 km/s/Mpc 
0.0432 
0.206 

-0.063 < nk,o < 0.017 
9.85 X 10^27 kg/in3 
4.25 X 10^28 kg/in3 
2.45 X 10^27 kg/n^3 



Table 2 A summary of result for scalar field power-law cosmology. Times are shown in Gyr for comprehensibility. Positive 
and negative fis's correspond to open and closed universes, respectively. For the interaction case, ao is set to 1. 





WMAP5+BAO+SN 




WMAP5 


a 




0.985 




1.01 




^k,Q = -0.0049 


-0.0179 < flkfi < 0.0081 


r^fc^o = -0.023 


-0.063 < flkfi < 0.017 


ao (non-int. case) 


1.9 X 10^^ 


ao > 9.85 X 10^6 (closed) 
ao > 1.5 X 10^'' (open) 


8.4 X 10^6 


ao > 5.1 X 10^*^ (closed) 
ao > 9.8 X 10^6 (open) 


non-interaction case 










^intercept 


2.65 Gyr 


2.62 Gyr < t < 2.69 Gyr 


2.65 Gyr 


2.55 Gyr < t < 2.76 Gyr 


^max 


3.99 Gyr 


3.94 Gyr < < < 4.01 Gyr 


3.97 Gyr 


3.82 Gyr < t < 4.13 Gyr 


^inflection 


5.33 Gyr 


5.27 Gyr < t < 5.40 Gyr 


5.29 Gyr 


5.09 Gyr < t < 5.50 Gyr 


interaction case 










^intercept 


2.77 Gyr 


2.73 Gyr < < < 2.80 Gyr 


2.76 Gyr 


2.66 Gyr < t < 2.87 Gyr 


^niax 


4.16 Gyr 


4.10 Gyr < t < 4.21 Gyr 


4.12 Gyr 


3.97 Gyr < t < 4.28 Gyr 


^inflection 


5.55 Gyr 


5.48 Gyr < t < 5.62 Gyr 


5.48 Gyr 


5.28 Gyr < t < 5.70 Gyr 
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potential and its first term in Fig. [T] eventually crosses, 
but it occurs much, much later at t — 8.8 x 10^^ Gyr- 
Along with the scalar potential, we also obtain 



Dc^ 



(14) 



in SI units. Using WMAP5+BA0+SN dataset and 
positive root, 



1.05 • 10^6 1.51 • 10 



23 



1.97 



3.01 • 10 



42 



1/2 I 

(15) 



where, for WMAP5 dataset alone, 



1.09 • 10^6 7.67 • 1024 9.37 • 10^^ 
+ - 



1/2 



^2 



^2.03 



^3.04 



dt.l 



In the late post-inflection phase, the first term is dom- 
inant over the k and D terms then the last two terms 
of the radicands are negligible (Fig. [I]) . The above two 
equations are approximated to find exact solutions. 



1/2 



hit 



where (po = [2aMpc/nj Into- Hence numerically 
using the WMAP5+BA0+SN combined dataselQ 

(j){t) w 4.17 X 10" - 1.03 X 10" Int 



whereas, for WMAP5 dataset alone. 



<t){t) « 4.23 X 10^ 



1.04 X lO'-'lni 



(18) 



(19) 



The radicand in ( 16 ) of the WMAP5 dataset is zero 



at approximately tinfloction — 5.3 Gyr, therefore so does 
(jiit). While the WMAP5+BA0-I-SN combined dataset 



has the zero radicand (hence zero 4'{t)) in (151 later at 



approximately t = 5.4 Gyr. With post-inflection ap- 
proximation, t w to exp {(j)/ y^2aMpc/ h) , the potential 
is written as function of the scalar field as 



y(0) 



3a2 



^2 
''0 



^-14,l^1aMlcln 



I ^fe'^^ ^-2v/^0/^/2Mgc/n 



^Scalar field exact solutions for the power-law cosmology with 
non-zero curvature and non-zero matter density are reported in 
l |Gumjudpai|2009^ 



Hence, we can plot V{4)) as 

V{ct)) « (2.72- 10-12) ei-92-io-"0 

+ (5.51- 10-1°) ei-9^-i°"'^ 
-(1.15- 10-1°) e2-88-io-^'^ 



(21) 



for WMAP5-hBA0-FSN dataset, whereas, for WMAP5 
dataset alone, 

V{ct)) « (5.94- 10-1°) ei-°2-io"''^ 



(1.36- 1024) e-i-9^ i°" 
(1.11 -1026) e-i-92.10- 



(22) 



The scalar equation of state as function of time can 
be found directly using (|3| in = p<^/p0, i.e. for 
WMAP5-hBA0-|-SN combined data 



(16) w<p = -0.333 + 



(4.43 • loi**) - imt 

(1.33 -1019) -(6.9 -102)^0. 



0.955 



955 _ ^0.985 



(23) 

lim wAt) = - 0.333, wJta) = -0.446, (24) 

f oo 

and for WMAP5 data alone 

(1.94-101^)- (6.25 • 10-4) ti oi 



(17) = -0.342 



(5.67 • 1017) - (6.97 • 10-2) ^l oi - t^M 

(25) 

lim wAt) = - 0.342 , wJto) = -0.452 . (26) 

The value of equation of state diverges when time sat- 
isfies condition. 



2 h Dctl"' kc'^tf' 



3M2 al t3"-2 al i2"-2 



(27) 



We find weighted effective value of equation of state 
parameter for dust and scalar field densities. That is 



PmC2 + P0C2 ' 

which diverges when time satisfies 

2 



t' 



(28) 



(29) 



Weighted value for WMAP5-HBA0-HSN combined data 
is 



(20) = -0.333 



'(1.87 • 104) -f7.01i' 



0.955 



(3.01 • 104) + (6.9 • 102) tO-955 + ^0.985 ' 

(30) 

«^off(io) =- 0.323, (31) 
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and for WMAP5 data afonc 

-(1.72 • 10^) + (6.25 • 10--*) i^ oi 



w^s = -0.342 



(2.52 • 103) + (6.97 • 10-2) ti oi + t^-o^ 

(32) 

wcff(io) =- 0.342. (33) 

These values do not match WMAP5 observation in 
wCDM model {w « —1). We will consider interaction 
case to check if effects from dark matter-dark-energy 
interaction could alter the equation of state and the 
scalar potential. 



5 Inclusion of interaction between dark energy 
and CDM 

In situation such that there is an interaction between 
the scalar field and cold dark matter which is a non- 
baryonic sector of dust ( Amendola||2000 1 , we can have 
interesting scenario. This can result from adding an 
Yukawa-like interacting term, —W {4>)mQ'il)ijj to the La- 
grangian density where tp is fermionic cold dark matter 



field as in (Das et. al. 2006) and (Amendola et. al 



2007). Other studies about dark sectors interaction can 



be seen in ( He et. al.|2007 ) The fluid equations for cold 
dark matter and the scalar field are then 



PCDM + 3i?pcDM 

+ 3Hp^{l + w^) 



TpcDM , 
— rpcDM : 



(34) 
(35) 



where the interaction rate is defined in term of the 
Hubble rate, F = 5(a)H or of the field speed F = Qcf) 
(( |Amendola||2000[ ) and ( [Wetterich ||1995[ )) and 5{a) is 
function of scale factor. This function relates to mass, 
of the fermionic CDM as 



5{a) 



d(ln m^{a)) 
d(lna) 



(36) 



or that is to say 'm^{a) — exp [J (5(a)d(lna)] . The mass 
function W{(j)) = m^{4){a))/mo in the Lagrangian is 
hence related to S and F as 



d{a) 



and 



d(lnVF(a)) 
d(lna) 



S{a)H 



1 dW ■ 



(37) 



(38) 



where toq is the present mass of the fermion field. 
Therefore Q = {l/W){dW /dcj)) and this can be inte- 
grated to obtain 



a /dW\ 

^^w[-d^ ' 



(39) 



and we found that 

s 



W{a) 



(40) 



assuming constant S. The interaction between dark 
matter and scalar field affects CDM matter density. It 



is straightforward from (34) to find that 



PCDM — PCDM.O I — 

V a 



3-5 



The equation ^ is also modified to 



■ Pcdmc 



1 dW\ dV 
W^) ~ d0 



(41) 



(42) 



After integrating the first term on right hand side, the 
effective scalar potential is found to be 



6-3 



PCDM.O C 



a 



3-S 



F(0). 



(43) 



The CDM density term is modified due to the interac- 
tion while the baryonic matter does not interact with 
the scalar field (this happens in chameleon scenario 
( Khoury~||2004 )). Hence in this situation, the equation 
( 14 ) is modified to 



kc^ 



(44) 



as well as the potential 



V = 



H 
"3 



3a2 



PCDM,0C 
203-5 

2a3 ■ 



(45) 



With post-inflection ap?proximation, the potential is 
written as function of the scalar field as, 



2kc^ 



3a2 



-20/^2ctM|c/;i 



PcBM.a c? 



2a: 



3-5 



/S(3-(5)0/^2M|c/S 



2al 



(46) 



Dark energy density that decays to dark matter density 
is Ap0 = pcDM,o/a^ - PcDM,o/a^"'' ■ For a positive 5, 
the scalar field decays into dark matter hence Ap^ < 
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while negative 6 gives reverse process. We define ef- 
fective scalar field density as p^f = P<p + ^Ptt> hence, 



„eff 



P0 = P<l> 



PCDM,0 



(1 - a^) 



(47) 



To find scalar field effective equation of state parameter, 
w'? , we consider effective scalar field fluid equation, 



.off 



= -iHpfil + wf) 



(48) 



Therefore, I + wf ^ (l/pf ) {P4> + Ap^ + + Ap^). 
Note that change in scalar field density, Ap^, results in 
change of scalar pressure density, Ap^. With the fact 
that, Ap0 = u'0(Ap0)c^, therefore indeed w'^^ = w^, 
i.e. unchanged. This is because of change in density 
would equivalently create change in pressure (this defers 



from (Das et. al. 2006j)). In order to account for total 
effective equation of state of all cosmic fluids, we define 
weighted equation of state parameter as 



Weff 



P4, + w^{Ap^)c 



(49) 



with the subscript eff instead of superscript. We pursue 
calculating all quantities as done in the non-interaction 
case. At present t — to, for interaction case we set oq = 
1 with constant interaction rate S — —0.03 obtained 



from (Guo et. al. 20071. Using WMAPS-^BAO-^SN 



combined dataset, the potential (as function of time) is 



V = 



1.03 • 10^^ 1.51 • 10^-^ 2.49 • 10*1 4.16 • lO'' 



^2 ' ^1.97 ^2.95 ^2.98 '| 

(50) 



whereas, for WMAP5 dataset alone, 
^_1.11-102'^ 7.67-1024 8.12-1042 1.33.1044 



2.03 



^3.04 



^3.07 



(51) 



in SI unit. The field solution for WMAP5+BA0+SN 
combined data is 



to 



1.06 - lO^*' 1.51 - 10^-^ 4.98 - 10 



-i41 



I ^2 ^1.97 ^2.95 

, \ 1/2 

8.33- 1042 \ . 

(52) 

(j){t) w 4.17 X 10" - 1.03 X 10"lni, (53) 



and for the WMAP5 alone, 




7.67-1024 1.63-1043 

£2.03 £3.04 
1/2 

dt, (54) 



2.66 - 10 



44 



^3.07 

(j){t) w 4.23 X 10" - 1.04 X 10^3 ^ (55) 
with estimation that the barotropic density and spatial 



curvature term are subdominant as before. Using (53) 



with WMAP5+BA0+SN dataset, the scalar potential 
is 



.2.72 - io-i2ei-92-io-"* 
-f 5.51-10-iV-95 io-"^ 
-1.91-10-"e2.88.io--^ 



9.61 - 10 



-llg2.90-10" 



(56) 



and using (55 1 with WMAP5 alone, the potential is 



y(</)) «5.94-10-i"ei-92.io— ^ 

+ 1.36-10-"ei-94i°" 
- 1.91 - 10-"e2.92.io- 

-9.12-10-"e2-95io" 



(57) 



The equation of state for the scalar field and weighted 
value with dark matter interaction as function of time 
can be found directly as previously. These are in form 
of complicated functions plotted in Figs. |2]and[3j The 
values at present are 



u'0(to) = -0.446, Wcff(io) 



-0.323, 



(58) 



for WMAP5+BA0+SN data and for the WMAP5 
dataset alone they are 



w^to) = -0.452 , Weff(io) = -0.342 . 



(59) 



For both datasets, w^{t) and Wcs{t) converge to ap- 
proximately — 1/3. In post-inflection period, non- 
interaction and interaction cases are less and less differ- 
ent due to the subdomination of the dust density term. 
These values at present time for non-interaction and in- 
teraction cases are the same at three significant decimal 
digits. The value of diverges when time satisfies the 
condition. 



fic I PCDM,0 ^0 



— cx5 



3a 



3M, 



Pbfi 



fcc2 ^2" 

al t2"-2 



(60) 
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Non-Interacting Case 
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WMAP5+BA0 + SN - 
WMAP5 ■ 




1x10^'' 2xl0" 3x10^^ 4xl0" 5x10^^ 6xl0" 
Interacting Case 

4 Gyr 8 Gyr 12 Gyr 16 Gyr 20 Gyr 



3. X 10"^ 




\\ WMAP5+BA0+SN 

\\ WMAP5 


2. X 10"^ 




a V(t) (inverse square term only) 


1. X 10"^ 











0{t) ' 




: 5.5 Gyr (post-inflection period 


1. xlO-5 








Ix 


10^'' 2x10^' 3x10^^ 4x10^' 5x10^^ 6xl0" 



Fig. 1 Potential plots for non-interaction case (upper 
panel) and interaction case (lower panel). The potentials 
are shown in solid lines. The units of the horizontal and 
vertical axes are sec (or Gyr) and J/m^, respectively. Dots 
mark their inflection points. The light-shaded regions cor- 
respond to ranges of Qk,o, and the central solid lines plot- 
ted using flk,o- Dashed lines are the potentials' first terms, 
showing that they are dominant at late time. Plotted in 
dotted lines are (p^{t) which is positive after tinflcction after- 
which <jf>(i) becomes real. In addition, each potential does 
not actually converge to its first term, but later intersect 
with and deviate from it, although they are still very close 
together (not shown). However, this occurs much later (at 
t = S.Sx 10®'' Gyr in both cases). The potential is lower in 
presence of the interaction between dark matter and dark 
energy. 



and for Wcs, it diverges when time satisfies condition, 
^2(a-i) _ —k{ctQ /aao)^ as in non-interaction case. 



Considering mass potential function ( 40 ) in scalar-field 
power-law cosmology, it is simply W = Wq (t/to)"^ . 
With similar procedure as used in finding V{(t>), we ob- 
tain 



Vr(0) « Wo exp 
^ \/a S 



(61) 
(62) 



The mass potential function is exponentially increasing 
if (5 > 0, i.e. when dark energy decays into dark matter. 
On the other hand, in our case, 5 = —0.03, we have 
exponentially decay mass potential function. The result 
for both datasets are 



Q«2.88x 10-l^ 

for WMAP5+BA0+SN data and 
l¥(0)«e2-92xio-^=0^ 
Q w 2.92 x 10-^^ 



(63) 
(64) 



(65) 
(66) 



using WMAP5 data alone. We set Wq = W{(j){ao)) = 1. 
W{(t)) is plotted in Fig. [i] 



6 Conclusion 

We consider scalar-field power-law cosmology featuring 
homogeneous canonical scalar field and dust fluid evolv- 
ing in a late non-flat FLRW universe expanding with 
power-law function. The scalar field is minimally cou- 
pled to gravity and dust fluids are baryonic and CDM 
perfect fluids. We use two sets of observational data, 
combined WMAP5-t-BA0-|-SN dataset and WMAP5 
dataset, as the inputs. Mean values of both sets suggest 
slightly closed geometry. With k = -t-1, it gives mean 
value of ao at 8.4 x 10^^ for WMAP5 and 1.9 x 10^^ 
for WMAP5-|-BA0-fSN datasets. For closed universe, 
the WMAP5 dataset puts the lower limit of 5.1 x 10^'^ 
for ao while the WMAP5-hBA0-hSN dataset dataset 
puts the lower limit of 9.85 x lO^^. The WMAP5 data 
alone yields the exponent a = 1.01 agreeing with the 
previous study of H{z) data which gives a = 1.07 (Dev 
et. al.||2008|). However combined WMAP5+BA0-hSN 



data disagrees with the previous combined i?(z)-fSN 
la data, i.e. WMAP5-I-BA0+SN gives a = 0.985 while 
H{z)-\-W gives a = 2.28 (closed geometry) ( |Dev et. ah 



2008^ 



Non- Interacting Case 



Hon- Interacting Case 

123456789 10 



^^^^^^^^ 5.3 Gyr {post-inflection 


lim iv^(t) 










ff WMAP5+BA0+SN 






WMAP5 



5.0x10-''' 1.0x10-" 1.5x10-" 2.0x10" 2.5x10-" 3.0x10-" 
Interacting Case 
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^^^^^^^^^^ 5.5 Gyr {post-inflection -*) 


lim iv^{t) 










'f WMAP5+BA0 + SN 






WMAP5 



5.3 Gyr {post-inflection -») 



lim Weff{t) {dashed) 



rveff { t ) { solid) 



WMAPS+BAO+SN - 
WMAP5 - 



5.0x10" 1.0x10^' 1.5xl0-^'' 2.0x10-^ 2.5x10^' 3.0x10^' 
Interacting Case 





1 5.5 Gyr {post-inflection -*) 


lim Weff{t) {dashed) 


^^eff (t) {solid) 




WMAP5+BA0+SN 

WMAP5 



5.0x10^^ 1.0x10^' 1.5x10^' 2.0x10^' 2.5x10^' 3.0x10 



5.0x10" l.Oxlo" 1.5x10" 2.0x10" 2.5xlo" 3.0xlo" 

Fig. 3 Weighted effective equation of state parameter, ujoff 
for non-interaction and interaction cases. Interaction effect 
lifts the values of Wefi at earlier time before entering post- 
inflection era when interaction effect is less important due 
to subdomination of the dust term. In both non-interaction 
and interaction cases, Wcff diverges when time satisfles con- 
dition, = -fc(cto/aao)^- Notice that here fc > 0. 
Existence of divergency in Wcfi, (condition that t is real) 
depends on the value of a — 1. 



Fig. 2 Equation of state parameter, for non-interaction 
and interaction cases. The dash lines are asymptotic value 
at late time. Horizontal axis are in sec and Gyr. The light- 
shaded regions correspond to the ranges of f2fc,o, and the 
central solid lines to A^.q. Interaction affects in slightly 
lowering the value due to decaying of dark matter to 
dark energy. In post-inflection period, their values (non- 
interaction and interaction cases) are less and less different 
due to the subdomination of the dust density term. 



Interacting Case 



2 10 
s 



5 



WMAP5+BA0+SN 




WMAP5 






during post-inflection period 


2x10" 


4x10" 6x10" 8x10" 1x10" 



* 

Fig. 4 Mass function VK(0) of the dark matter-dark en- 
ergy interaction term versus 0. The mass function declines 
sharply at flrst and slowly converges to zero, f inflection cor- 
responds to </> 4 X 10^* for both datasets. 
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The scalar field potential is found with the two 
datasets. Slope of the potential is inflected after 
^inflection- This time scale characterizes the time when 
the Hubble expansion and H terms begin to dom- 
inate over the dust and curvature terms. For the 
WMAP5 dataset, iinflcction =5.29 Gyr and for the 
WMAP5+BA0+SN dataset, tinflcction = 5.33 Gyr. We 
also consider situation when there is interaction be- 
tween dark energy-dark matter in order to see its ef- 
fects on the potential and equation of state. Con- 
stant interaction rate is assumed here with 5 = —0.03 



from (Guo et. al. 2007) corresponding to dark mat- 



ter decaying into the scalar field. For interaction case, 
^inflection =5.48 Gyr (WMAP5 dataset) and ^inflection = 
5.55 Gyr (WMAP5-hBA0-hSN dataset). The interac- 
tion affects in lowering the height of scalar potential and 
shifting the potential rightwards to later time (see Fig. 
[l]and Table[2|. In this study, with closed geometry sug- 
gested by WMAP5 data, the field equation of state does 
not match the WMAP5 observation (wCDM model), 
i.e. w « — 1. Effects of interaction is in slightly lowering 
the Wtp value while slightly lifting up the weighted effec- 
tive equation of state Wcs at earlier time. We found that 
the mass potential function W{4>) of the interaction is 
approximately in form of exponential decay function. 
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